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Introduction 


Structural and energetics properties of atoms near a grain 
boundary are of great importance from theoretical and experimental 
standpoints. From various experimental work it is concluded that 
diffusion at low temperatures at polycrystalline materials take 
place near grain boundary. Experimental and theoretical results 
also indicate changes of up to 70% in physical properties near a 
grain boundary. The Embedded Atom Method (EAM) (1) calculations on 
structural properties of Au twist grain boundaries (2) are in quite 
good agreement with their experimental counterparts. The EAM is 
believed to predict reliable values for the single vacancy 
formation energy as well as migration energy. However, it is not 
clear whether the EAM functions which are fitted to the bulk 
properties of a perfect crystalline solid can produce reliable 
results on grain boundaries. 

One of the objectives of this work is to construct the EAM 
functions for Cu and use them in conjuction with the molecular 
static simulation to study structures and energetics of atoms near 
twist grain boundaries E5 and E13 in Cu. This provide test of the 
EAM functions near a grain boundary. In particular, we determine 
structure, single vacancy formation energy, migration energy, 
single vacancy activation energy, and interlayer spacing as a 
function of distance from grain boundary. Our results are compared 
with the available experimental and theoretical results from grain 
boundaries and bulk. 


Methodology 

Interatomic pair potentials suffer from two major problems, 
zero Cauchy presure, i.e., C 11 -C 12 =0, and equality of single vacancy 
formation energy E 1v and cohesive energy E c . In practice, for most 
crystals and E 1v *E c .To overcome these and other shortcomings 
of pair potentials the EAM potentials are developed. In the EAM, 
the major contribution to energy of each atom is due to the 
embedding energy which is supplemented by a two body potential . The 
embedding energy of an atom can be interpretted as the energy that 
is required to embed that atom in the electronic charge produced by 
the other atoms in the crystal . The two body potential <p can also 
be interpretted as the electrostatic interaction between cores of 
embedding atom and atoms of crystal. Energy of atom i in the EAM 
formalism can then be written as, 


E i = ^(pj) +0.5 


£.<t> 




i i ^ ij ^ 


where E. is the energy of atom i, F 1 is the embedding energy of atom 
i, <p Vi i's the two body potential between atoms i and j, r^- is the 
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interatomic distace between atoms i and j, and is the charge 
density at site i due to atoms in other sites. One of the 
assumptions in the EAM approach is that p i is superposition of 
atomic charge densities from all the other sites. 


Pi 


• Z. 


P* 


where p a is the atomic charge density. 

F and <p are determined by considering for them functional 
forms and finding their parameters by fitting to the bulk 
properties of crystalline solid (3) . 

In our calculations, we have employed the EAM functions of Cu 
along with the molecular static routine to study structures and 
energetics of twist grain boundaries. In a molecular static 
simulation, total energy of the system is minimized with respect to 
positions of atoms. In any energy minimization approach extreme 
caution should be taken that final state of the system has the 
lowest energy (i.e., global minimum). 


Results 

a) single vacancy formation energy 

Our 25(6=36.87) and 213(0=22.62) twist grain boundaries are 
comprised of 14 fee (001) atomic layers, 80 atoms per layer, and with 
half of the layers rotated with respect to the (001) axis by an 
angle teta. A vacancy is created in a desired site and the single 
vacancy formation energy Elvf is calculatedusing. 


Ef v = E(N- 1,1) - E(N, 0) 


where E(N-1,1) is the total energy of a system of N atomic sites 
with N-l sites filled and 1 site vacant and E(N.O) is the energy 
with all N sites filled. Our results for E. * are plotted in 
Fig. (2). 


b) Migration and activation energies 

Various migration paths for vacancy is considered (look at Fig. 
1) by moving the atom through the desired path and allowing it to 
relax in the plane perpendicular to the path. All the other atoms 
are allowed to relax without any constraint. Total potential energy 
of the system is then calculated as a function of migrating 
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distanse. Our results for migration energies are reported 
below and plotted in Fig. (2): 
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in Fig. (2) . 

c) Interlayer spacing across a 25 grain boundary 

Total energy of the grain boundary lattice described in part (a) 
is minimized using our molecular static code and interlayer spacing 
is calculated and plotted in Fig. (3). 


Summary and conclusion 

Our work can be summarized as follows: 

a) Vacancy formation energy is shown to be directly correlated to 
the degree of coincidence. 

b) Single vacancy activation energy near the grain boundary is 
about 60% lower than its corresponding value in the bulk. 

c) Vacancy migration energy near a 25 or 213 grain boundaries is 
about 0.1 ev which is much smaller than its corresponding bulk 
value of *0.7 ev. 

d) Interlayer spacing at a 25 grain boundary is enhanced by a 
factor of about 17% as compared to its bulk value. 
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